Dynamics of the breakdown of granular clusters 
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Recently van der Meer et al. studied the breakdown of a granular cluster (Phys. Rev. Lett. 88, 
174302 (2002)). We reexamine this problem using an urn model, which takes into account fluctua- 
tions and finite-size effects. General arguments are given for the absence of a continuous transition 
when the number of urns (compartments) is greater than two. Monte Carlo simulations show that 
the lifetime of a cluster r diverges at the limits of stability as r ~ N 1 * 13 , where N is the number of 
balls. After the breakdown, depending on the dynamical rules of our urn model, either normal or 
anomalous diffusion of the cluster takes place. 



I. INTRODUCTION 

Dissipation of kinetic energy during inelastic colli- 
sions in gaseous granular systems has profound conse- 
quences jl], |^] . One of the most spectacular ones is forma- 
tion of spatial inhomogeneities || , which drastically con- 
trast with a uniform distribution of molecules or atoms 
whose dynamics is essentially elastic. 

Some time ago Schlichting and Nordmeier presented 
a simple experiment which demonstrates some conse- 
quences of inelasticity of granular systems 0). They 
used a container separated into two equal compartments 
by a wall which has a narrow horizontal slit at a cer- 
tain height. The container is filled with balls (plastic 
or metallic) and subjected to vertical shaking. For vig- 
orous shaking the balls distribute equally between two 
compartments. However, when the shaking is sufficiently 
mild, a nonsymmetric distribution occurs. In such a case 
the compartment with majority of balls, due to numer- 
ous inelastic collisions, is effectively cooler than the other 
one. Consequently, less balls are leaving this compart- 
ment which stabilizes such an asymmetric distribution of 
balls. To explain this experiment, Eggers derived a phe- 
nomenological equation for the flux F(n) of balls leaving 
a given compartment || 



F(n) = Cn 2 exp{-Bn 2 ). 



(1) 



In the above equation n is the concentration of balls in 
a given urn and B and C are constants which depend 
on the properties of balls, typical sizes of the system and 
of parameters of shaking (the constant C may be elim- 
inated by an appropriate redefinition of the time scale) . 
In agreement with experiment, eq. (|l]) predicts for suffi- 
ciently large B unequal distribution of balls. The above 
experiment was repeated in the case when the number of 
compartments L was greater than two by van der Meer 
et al. |(| . In such a case formation of unequal distribu- 
tion of balls is accompanied by strong hysteresis which 
is in agreement with theoretical analysis Q. Moreover, 
certain aspects of these phenomena for L = 2 were ap- 
proached using hydrodynamic equations Q . 

Recently, van der Meer et al. examined the case of 
L > 2 further 0] . In particular, they studied dynamics of 
configurations (clusters) starting from all balls localized 



in a single compartment. Using a theoretical model based 
on eq. (Q), they have shown that when shaking is strong 
enough such a cluster breaks down and diffuses with the 
anomalous diffusion exponent 1/3 (in the following we 
refer to this model as MWL). For less vigorous shaking, 
the cluster remains relatively stable and only after some 
time it abruptly breaks down. Some of their predictions 
were confirmed experimentally. 



In the framework of the MWL model it is rather dif- 
ficult to include the effect of fluctuations. Such fluctua- 
tions might originate due to for example a finite number 
of balls and especially close to critical points they might 
play an important role. In an attempt to take such ef- 
fects into account a generalization of Ehrenfest' s |l0[ urn 
model was recently examined in the case L = 2 Jll|. Rel- 
ative simplicity of the model allows for a detailed study 
of its various characteristics. 



The motivation of the present paper is to re-examine 
the breakdown of granular clusters using the urn model 
in the case L > 2. In section [il] we define the model 
and present its steady-state phase diagram for L = 3. 
We also argue that, in analogy to the Potts model in 
the mean- field limit, there are no continuous transitions 
for L > 2. In section [II we examine dynamics of the 



breakdown of clusters in a similar way as van der Meer 
et al. H ■ Although qualitatively our results are similar to 
theirs, in our model the diffusion of the cluster is normal 
with the exponent 1/2. Moreover, we calculate the size 
dependence of the lifetime of a cluster r and show that 
at the limits of stability it scales as iV 1 / 3 . In section IV 
we present a modified version of the urn model which in 
the steady state reproduces the flux ([[]). The diffusion 
of the broken down cluster is then shown to be anoma- 
lous with exponent 1/3, as it was already found ||. It 
was suggested that essential features of the MWL model 
are independent on the precise form of the flux (jl|), as 
long as it has a single hump |)|. On the contrary, our 
results show that at least the diffusion exponent depends 
on some details of the flux and not only on its qualita- 
tive shape (in our models the flux is also a single hump 
function). Section |v| contains our conclusions. 
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II. MODEL AND ITS STEADY-STATE 
PROPERTIES 

Our model is a straightforward generalization of the 
two- urn case jll| : N particles are distributed between L 
urns and the number of particles in i-th urn is denoted 
as Ni (J2i=i W = N). Urns are connected through slits 
sequentially: i-th urn is connected with [i — l)-th and 
(i + l)-th. Moreover, periodic boundary conditions are 
used, i.e., first and L-th urns are connected. Particles in 
a given urn (say i-th) are subject to thermal fluctuations 
and the temperature T of this urn depends on the number 
of particles in it as: 



0.25 



T(m) = T + A(l - m), 



(2) 



where Tij is a fraction of the total number of particles 
in a given urn (rij = Ni/N) and T and A are positive 
constants. Equation (0) is the simplest function which 
reproduces the fact that due to inelastic collisions be- 
tween particles, their effective temperature decreases as 
their number in a given urn increases. Next, we define 
the dynamics of the model as: 

(i) One of the N particles is selected randomly. 

(ii) With probability exp[— 1/T(rii)] the selected parti- 
cle is placed in a randomly chosen neighboring urn, 
where i is the urn of a selected particle. 

The above rules implies that the flux of particles leav- 
ing i-th urn is, up to a proportionality constant, given 
by 



F{m) = rn exp 



T(r, 



(3) 



where T(ni) is defined in (^). Let us notice that the flux 
(||), similarly to ([!]), is a single hump function. Having 
an expression for the flux we can write the equations of 
motion as: 

_l = l F (n i _ 1 ) + i F {n i+1 )-F{n i ), (4) 

where i — 1,2,..., L. Steady-state properties of this 
model can be obtained using similar analysis as in the 
L = 2 case or as for L > 2 but with fluxes given by 
cq. ([!]) . The results of this analysis in the L = 3 case 
are presented in Fig. [l]. In region I and II the symmetric 
phase (ni =712 = 713 = 1/3) is stable. The continuous 
line in Fig. [ij which locates the limit of stability of this 
phase, is given by the following equation: 



_ 2A 
Tn = \ — — . 



(5) 



This equation has a very similar form to the correspond- 
ing equation in the L = 2 case |Tj]| . Asymmetric so- 
lution, where one of the urns has the majority of balls 
and remaining two urns have only a small equal frac- 
tion of balls (n 1 > ri2 = n^), is stable in region II and 
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FIG. 1: The steady-state phase diagram for the three-urn 
model. See text for a description of phases. 



III. The line separating regions I and II can be deter- 
mined only numerically as a solution of a transcendental 
equation, similarly to the L = 2 case |ll[. There is also a 
third type of solution where two urns contain majority of 
balls and the third urn has only a small fraction of them 
(rii = ri2 > 713). Such a solution, which has saddle-like 
stability, exists only in region III. Similar solutions can be 
found for the MWL model 0, fjj. An important, qualita- 
tive difference with the case L — 2, is that regions I and 
III are always separated by region II where both sym- 
metric and asymmetric solutions are stable, hence the 
tricritical point is located at the origin To = A = 0. It 
means that a phase transition between these two phases 
is always accompanied by hysteresis effects. On the other 
hand in the L = 2 case continuous transitions are possi- 
ble, which are not accompanied by hysteresis JIl[ . Such 
a behaviour is actually in agreement with experimental 
data and with MWL model ||. 

Has this qualitative difference a more general expla- 
nation or is it rather a coincidental property? In our 
opinion, absence of continuous transitions for L > 2 is 
a generic property of such systems and at least to some 
extent could be understood. First, let us notice that 
the phase transition for L = 2 is a manifestation of the 
spontaneous symmetry breaking in the system: in certain 
regime one of the two identical urns is preferentially filled 
with balls. Such a situation resembles the phase transi- 
tion in the S = 1/2 Ising model, where below certain 
temperature the up-down symmetry is broken and the 
system acquires spontaneous magnetization Actu- 
ally, this analogy can be confirmed more quantitatively. 
We have shown that for L = 2 and at the critical point 
the probability distributions has the same moment ratios 
as in the Ising model in dimension d greater than the so- 
called upper critical dimension (d > 4) 113]. Let us no- 
tice, that in our model balls are selected randomly which 
means that this is essentially a mean-field model. More- 
over, our model is a dynamical, spaceless model, contrary 
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to the Ising model, which is a lattice equilibrium model. 
The fact that such different models have some similarities 
shows that as far as the critical behaviour is concerned 
what really matters is symmetry. In both cases this is the 
Zi symmetry which is broken below the critical point. 

Pushing this analogy further, we expect that for L > 2 
the phase transition in our model should be similar to 
the phase transition of the L-state Potts model above 
the critical dimension juj. In the L-state Potts model 
at sufficiently low temperature one of the L symmetric 
ground states is preferentially selected. However, it is 
well-known that above the upper critical dimension and 
for L > 2 there are only discontinuous transitions in the 
Potts model |Q . Consequently, the transition in the urn 
model, and most likely in related models which preserves 
Zl symmetry of compartments, should be discontinuous. 

Let us notice that one can easily break the symmetry 
of the compartments e.g., changing the boundary condi- 
tions, which in our analogy introduces some asymmetry 
in the Potts model. It is possible that in such a case 
the system effectively will become similar to the L = 2 
system and will exhibit a continuous transition. Finally, 
we expect that for L > 3 the phase diagram should be 
topologically similar to the one for L = 3 shown in Fig. [I]. 



III. DYNAMICAL PROPERTIES OF CLUSTER 
CONFIGURATIONS 

In the present section we study certain dynamical 
properties of cluster configurations. We used Monte 
Carlo simulation. Since it is rather straightforward, we 
omit a more detailed description of the numerical imple- 
mentation of the dynamical rules of our model. Initially, 
we place all balls in one urn and examine its subsequent 
evolution. If the parameters To and A are such that the 
system is in region I then such a cluster is unstable and 
after some time due to fluctuations it breaks down and 
balls spread throughout all urns. This is illustrated in 
Fig. ^| which shows the concentration of balls in the urn 
in which the balls were initially placed. Let us notice 
that (i) the breakdown is relatively abrupt and during 
the evolution up to the breakdown the concentration of 
balls only slightly decreases; (ii) upon approaching the 
line separating regions I and II the lifetime of the cluster 
r increases. Since in region II the asymmetric state has 
an infinite lifetime it means that r must diverge upon ap- 
proaching this region. This behaviour is seen in Fig. |3j. 
In addition to the three-urn case we also made analo- 
gous measurements of r for L = 5 and 7 and the results 
are also shown in Fig. ||. Let us notice that results pre- 
sented in Fig. U and Fig. || are similar to those obtained 
by van der Meer ||, although they are parametrized by 
a different variable. 

The limit of stability of the asymmetric phase can be 
regarded as a critical point. Thus, we expect that exactly 
at this point e.g., the lifetime r has a power-law diver- 
gence r = N z , and z > 0. Such a behaviour is shown 
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FIG. 2: The time evolution of the fraction of balls of the 
cluster n c i close to the limits of stability of the asymmetric 
phase (N = 5 ■ 10 4 , L = 3). The values of T are indicated. 
For A = 0.3 the limit of stability of the asymmetric phase 
is at T = 0.169829772... For a larger number of balls N, 
stochastic fluctuations will diminish. 
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FIG. 3: The average lifetime of a cluster r as a function of 
To for different number of urns L. Each point is an average 
of at least 300 runs. 



in Fig. ^. From the slope of the straight line, which is a 
least-square fit to our data we estimate z = 0.32(3). Let 
us notice that in the two-urn model at the limits of stabil- 
ity t exhibits a very similar divergence |ll| . In the case 
L = 2 more precise calculations were possible strongly 
suggesting that 0=1/3 which is also consistent with the 
present three- urn model result. Let us emphasize that 
because in our model the number of balls is finite, we 
can study size dependent quantities as shown in Fig. |^. 
Such calculations would not be possible for models solely 
based on steady-state equations. 

Finally, let us examine the breakdown of a cluster in 
the many-urn case L » 1. In such a case a continuous 
approach to the MWL model shows that after breaking 
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IV. THE PAIR MODEL 

One can easily construct urn models for which the ex- 
pression for the flux will have a different form. In particu- 
lar, redefining the effective temperature (0) and drawing 
each time a pair of balls we obtain an urn model with the 
flux of exactly the same form as eq. (|f]). This dynamics 
takes into account some of the two particles correlations. 
It allows us to recover some properties of the MWL model 
and establish further results. 

The model, which we call a pair model, is similar to 
the previously described one, except that its dynamics is 
now defined as: 

(i) Two different balls are selected randomly. 



FIG. 4: The average lifetime of a cluster r as a function of the 
number of balls N at the limits of stability of the asymmetric 
phase. Each point is an average of at least 300 runs. 
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FIG. 5: The average occupancy of a central urn N c i as a 
function of time t. The slope of decay is very close to 0.5 which 
confirms the diffusive nature of spreading (N c i ~ t -1 / 2 ). Each 
curve is obtained from averaging over 50 independent runs. 



down, the cluster diffuses with the anomalous exponent 
1/3 g. Results of our simulations are shown in Fig. 0. 
From these data we conclude that spreading of a clus- 
ter occurs with the ordinary exponent 1/2 rather than 
anomalously. Ordinary diffusion in our model can be 
also easily explained analytically applying basically the 
same continuous approach as used in H. In this approach 
the set of equations of motion (Q) is transformed into a 
partial differential equation. Then, one immediately re- 
alizes that the linear term in front of the exponent in 
eq. (||) leads to the ordinary diffusion equation. On the 
other hand, the anomalous diffusion of MWL model can 
be traced back to the quadratic (in n) term in the flux in 
eq. (|f]). This quadratic term is related with two-particle 
collisions ||. 



(ii) If and only if the two balls are in the same urn, 
with probability exp[— Bnf] the selected balls are 
placed in the same randomly chosen neighboring 
urn, where i is the urn of the selected particles. 

One can easily see that the probability that two ran- 
domly selected balls belong to the i-th urn is given as 
t£ • jfzzr i which for N — > oo becomes nf . Multiplying nf 
with the transition probability exp[— Bnf ] we obtain that 
the flux in the pair model is proportional to eq. (|f|). ft 
means that as far as the steady-state properties are con- 
cerned, the pair model is equivalent to the MWL |(| [?|] . In 
particular for L = 2 one easily obtains the critical value 
B = 4 for the continuous transition between the symmet- 
ric (B < 4) and asymmetric phase (B > 4). For L = 3 
one obtains two critical points Bi = 6.55270341 f . . . and 
£?2 = 9. The first one can only be determined numer- 
ically. Similarly to Fig. EL for B < B2 the symmetric 
solution is stable whereas for B > B\ the asymmetric 
solution stable. In the interval B S [B±,B2] both sym- 
metric and asymmetric solutions are stable, which is the 
interval showing hysteresis with respect to the driving 
parameter B. 

Qualitatively the dynamical properties of cluster con- 
figurations in the pair model are similar to those de- 
scribed in previous section. In particular for L = 3 and 
B = Bi, the average lifetime of a cluster r as a function 
of the number of balls N once more shows a power-law 
divergence r = N z , with z = 0.31(3) suggesting that 
z = 1/3. It shows a certain universality of this exponent 
with respect to different dynamical rules. 

Finally, Fig. |^ shows the diffusion of the broken down 
cluster. Since the asymptotic slope of our data is very 
close to 1 /3 we conclude that in this case the diffusion is 
anomalous, as already predicted by van der Meer et al. 
who used the continuous approach [§j . 

The pair model and the model examined in the pre- 
vious section exhibit qualitatively similar behaviour for 
most of the physical quantities. The main difference is 
the diffusion: it is anomalous in the pair model and or- 
dinary in model examined in the previous section. It 
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qualitatively recover experimental findings and previous 
steady-state calculations. In addition, our models take 
into account fluctuations caused by the finite number of 
balls. Using symmetry properties, we related them with 
high-dimensional Potts model and argued that for L > 2 
phase transitions in such systems should be discontin- 
uous. Although several quantities exhibit qualitatively 
a similar behaviour for the two different versions of the 
model, there are important differences too. In particu- 
lar, these models predict a different diffusion of a broken- 
down cluster, which could be cither ordinary or anoma- 
lous. It shows that the type of diffusion is very sensitive 
to dynamical rules of the model, and consequently, to the 
form of the flux. 



FIG. 6: The average occupancy of a central urn N c i as a 
function of time t for the pair model. The slope of decay 
is very close to 1/3 which confirms the anomalous diffusive 
nature of spreading (N c i ~ t -1//3 ). Each curve is obtained 
from averaging over 50 independent runs. 



would be desirable to experimentally examine the nature 
of diffusion in such systems. 
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